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Abstract. We present theoretical continuum emission spectra 
(SED's), isophotal maps and line profiles for several models 
of photoevaporating disks at different orientations with respect 
to the observer. The hydrodynamic evolution of these models 
has been the topic of the two previous papers of this series. 
We discuss in detail the numerical scheme used for these di- 
agnostic radiation transfer calculations. Our results are quali- 
tatively compared to observed UCHll's. Our conclusion is that 
the high fraction of "unresolved" UCHll's from t he cat alogues 
of Wood & Churchwell (|1989|) and Kurtz et al. (|1994 cannot 



be explained by disks around massive stars. In particular, the 
observed infrared spectra of these objects indicate dust tem- 
peratures which are about one order of magnitude lower than 
expected. We suggest that disks around close companions to 
OB stars may be necessary to resolve this inconsistency. Alter- 
natively, strong stellar winds and radiative acceleration could 
remove disk material from the immediate vicinity of luminous 
O stars, whereas for the lower luminosity sources considered 
here this will not occur. We also find that line profiles tracing 
the evaporated material originating from the disk are not influ- 
enced significantly by the existence of stellar winds over a wide 
range of wind velocities (400 - 1000 km s^^). We compare our 
results to the bright IRAS source MWC 349 A. Many of its 
properties, especially its spatial appearance in high-resolution 
radio maps, can be well explained by a disk surrounding a UV 
luminous star with a high velocity stellar wind. 

Key words: Radiation transfer - Line: profiles - Stars: circum- 
stellar matter - ISM: H II regions - ISM: jets and outflows 



1. Introduction 

A large fraction of early-type main sequence stars are as- 
sociated with ultracompact Hll-Regions (UCHlls, Wood & 
Churchwell 1989 1. These are characterized by large electron 
densities iVo > lO"* cm^"^, sizes < 0.1 pc and temperatures 



Te « 10 000 K. The overpressure in these regions should lead 
to expansion and dissipation on time scales of a few thousand 
years. Considering the expected lifetimes of massive stars of 
several million years, however, the high abundance of observed 
UCHlls translates into UCHll mean lifetimes of several 10^ 
years. This contradiction can be resolved by: 1) the UCHlls 
could be constrained by high pressure in thei r vicin ity, or 2) by 
gravitationally infalling material (Reid et al. 1981 ), or 3) there 
could exist a process which continuously "feeds" the UCHlls 
with matter High pressures certainly can be expected in the 
highly turbulent molecular cloud cores, which are the birth- 
places of young massive stars (De Pree et al. 1995, Garcfa- 
Segura & Franco |1996| , Xie et al. |1996| ). Still, it is not clear 
how turbulence in a cold clumpy medium can contain the warm 
(T ~ 10"* K), high density ionized material for extended peri- 
ods of time - many of the technical details of this proposal need 
to be worked out. 

The photoevaporating disk model proposed by HoUenbach 
et al. ( |1993| ) and Yorke & Welz ( |1993| ) offers an attractive al- 
ternative. A circumstellar disk around a luminous OB star is 
continuously photoionized by the central source. The existence 
of a powerful stellar wind can modify the quantitative details 
of this model, but the basic result remains the same. Long- 
lived UCHlls are the necessary consequence of disks around 
hydrogen-ionizing sources. In a subsequent paper by HoUen- 



bach et al. (1994) the quasi-steady state structure of disks a- 



round ionizing sources with winds has been calculated (semi-) 



analytically and in Yorke (1995), Yorke & Welz (1996, here 



after Paper I), and Richling & Yorke (1997, hereafter Paper II) 
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the evolution of such circumstellar disks has been followed nu- 
merically under a variety of conditions. 

In Paper I it has been stressed that the phenomenon of 
disks in the process of photoionization is not restricted to the 
(presumably highly symmetrical) case of circumstellar disks 
around OB stars. Disk formation is a common by-product of 
the star formation process. Because OB stars seldom form in 
isolation, close companions with disks to a powerful source of 
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ionizing UV radiation and a stellar wind should be common. 
Strongly asymmetric UCHlls should result. 



Wood & Churchwell (1989) observed 75 UCHlls at A = 



2 cm and 6 cm with spatial resolution of 0. 4 using the VLA 
telescope and classified them by their spatial morphological 
structure into several types: 

- cometary shaped (20%), 

- core-halo (16%), 

- shell type (4%), 

- iiTegular or multiply peaked (17%) and 

- spherical or unresolved (43%). 

In order to interpret these observations in light of the pho- 
toionized disk models, further work must be done in refining 
the hydrodynamical models for the asymmetric morphologi- 
cal configurations expected when a disk is ionized by external 
sources. Diagnostic radiation transfer calculations of these nu- 
merical models are necessary for a quantitative comparison. 

Goal of the present investigation is to determine spectral 
characteristics and to calculate the expected isophote maps 
of the symmetrical UCHlls which result from circumstellar 
disks around OB stars. We are restricted by the limited num- 
ber of star/disk configurations which have been considered to 
date. We discuss in detail the physical (Sect. |^) and numerical 
(Sect. |) models of radiation transfer which we employed. The 
results for selected hydrodynamical models from Papers I and 
II are discussed in Sect. ^ and compared to observations of spe- 
cific sources in Sect. || We summarize our main conclusions in 
Sect. |. 

2. The physical model 

In Papers I and II of this series the time dependent photoevapo- 
ration of a 1 .6 M0 circumstellar disk around a 8.4 M0 star was 
calculated under a variety of physical conditions. The ionizing 
flux of the central source and its "hardness" as well as the stel- 
lar wind parameters (mass loss rate and terminal velocity) were 
varied. States of these models at selected evolutionary times are 
the basis for our diagnostic radiation transfer calculations. 

2.1. Continuum transport 

To determine the continuum spectral energy distribution (SED) 
over a frequency range from the radio region up to the optical, 
we take into account three major radiation processes: thermal 
free-free radiation (i.e. bremsstrahlung of electrons moving in 
the potential of protons in the H ll-region), thermal dust radia- 
tion and the radiation emitted from the photosphere of an em- 
bedded source. 

2.1.1. Free-free radiation 

For this process w e adop t the approximation for the emission 
coefficient (Spitzer 1969|): 
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Table 1. Parameters for the grain species used in the dust model 
ofPreibisch et al. (|1993|). 



Grain Species 


r,ub/[K] 


log ad /[cm] 


lognd/[g ^] 


aC 


2000 


-6.024 


14.387 


Silicate 


1500 


-5.281 




CMP 


125 


-5.222 


12.218 



Here, A'c and iVp are the particle densities of electrons and pro- 
tons. All other symbols have their usual meanings. We approx- 
imate the Gaunt factor gg for a non-relativistic plasma by: 



ge 
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where 7 (« 0.577) is Euler's constant. Assuming the validity 
of Kirchhoff's Law Sy — e^/K^ — B^iTf,), the absorption 
coefficient for thermal free-free radiation can be written (hv ^ 

fcTe): 



4(27r)i/2e6Ar^7V, 



p5ff 



(3m/c)3/2cTe/2^2 



(3) 



2.1.2. Dust emission 

We adopt the 'dirty ice' dust model developed by Preibisch et 
al. (1993), which includes two refractory components: amor- 
phous carbon grains (aC) and silicate grains as well as volatile 
ice coatings on the surface of the silicate grains at temperatures 
below 125 K (Core Mantle Particles, CMP's). The icy coatings 
contain 7% of the available amorphous carbon and consist of 
water and ammonium with a volume ratio of 3:1. At tempera- 
tures above 125 K the silicate core and approximately 1 1 amor- 
phous carbon particles are released into the dusty gas for each 
CMP. In Table |l] the sublimation temperature Tgub, the mean 
radius Sd and the number of grains per gram gas rid are listed 
for the different species. 

The absorption coefficient [cm~^] for the individual dust 
components is given by: 



(4) 



where the mean absorption efficiency Q^j^ for grain type "d" 
has been determined using Mie theory for spherical grains of 
a given size distribution. Figure |l] displays the absorption ef- 
ficiencies for the different dust components as a function of 
frequency. Each dust component's contribution to the source 
function due to thermal emission is also calculated under 
the assumption that = Bi,{T^). 
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Fig. 1. Mean absorption efficiencies for the different dust com- 
ponents. Solid line: amorphous carbon, dotted line: silicate, 
dashed line: CMP's. 

2.1.3. Net continuum absorption and emission 

Both emission processes mentioned above occur simultaneous- 
ly within the same volume. Thus the net absorption coefficient 
and source function are: 



1 



(5) 



(6) 



2.2. Forbidden lines 

In order to calculate profiles of the forbidden lines for the ele- 
ments oxygen and nitrogen ([O ll] 3726, [O III] 5007 and [N ll] 
6584), we adopt the following procedure. First, the equilib- 
rium ionization structure of these elements is calculated over 
the volume of consideration. Next, the occupation densities of 
metastable levels A'i due to collisional excitation by electrons 
is determined. We take into account Doppler shifts due to bulk 
gas motions and thermal Doppler broadening to calculate the 
profile function </>: 
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exp 



~x 2 



with the thermal Doppler width: 



(7) 
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Here R is the gas constant, /i the atomic weight of the rele- 
vant ion, and i> = i/o(l + wr/c) is the transition frequency 
Doppler-shifted by the radial velocity vr of the gas relative to 
the observer 

The emission coefficient of the transition k — > j, which 
enters into the equation of radiative transfer, is then given by: 



47r 



(9) 



where yl^j is the Einstein coefficient for spontaneous emission. 
Note that we have neglected radiative excitation and stimulated 
emission in this approximation. 

2.2.1. Ionization equilibrium 

The equations for ionization equilibrium for two neighboring 
ionization stages, r and r + 1, are: 



(10) 



We solve these equations simultaneously for the iV' up to the 
ionization stage ?■ = 3 for both oxygen and nitrogen. 



Radiative ionization. The rate of radiative ionization is calcu- 
lated from the flux of incident photons and the absorption 
cross section al, integrated over all ionizing frequencies. We 
use the radiation field of the central source and neglect scatter- 
ing to determine f^: 



1 B,{T,)Rl^ 
= — exp(-T). 



(11) 
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An analytical expression for the absorption cross section (j'l is 



given in Henry ( 197C ) 



Collisional ionization. This ionization process is important in 
hot plasmas, where the mean kinetic energy of the electrons 
is comparable to the ionization potentials of the ions. N I, for 
example, has an ionization potential 14.5 eV; the correspond- 
ing Boltzmann temperature is ^ 170000 K. The coefficient for 
collisional ionization is approximated by the analytical ex- 
pression in Shull & van Steenberg (1982). 



Radiative recombination. This is the inverse process to radia- 
tive ionization. For the recombination coefficient Q!r we use 



the formula given in Aldrovandi & Pequinot (1973, 1976). 



Dielectronic recombination. The probability for recombina- 
tion is enhanced when the electron being captured has a kinetic 
energy equal to the energy necessary to excite a second elec- 
tron in the shell of the capturing ion. The density of excited 
levels in the term scheme of the ions grows with energy. Thus, 
this process becomes more and more important with increasing 
temperature. We use two analytical expressions for a-o: one 
for temperatures between 2000 K and 60 000 K (Nussbaumer 
& Storey 1983 1 and one for higher temperatures (Shull & van 
Steenbergll983). 



Charge exchange. The exchange of electrons during encoun- 
ters of atoms and ions, e.g. + + iJ+ is also 
important. Arnaud & Rothenflug ( 1985| ) give an expression for 
the coefficients J"'. Special care is necessary in the case of the 
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reaction O'^+H'^. Due to the similarity of the ion- 

ization energies of hydrogen and oxygen {AE = 0.19 eV) the 
backward reaction is also very effective. At sufficiently high 
electron temperatures this leads to the establishment of an ion- 
ization ratio Noo/Nq+ k, {9/8)Nifo /N[j+, even in the ab- 
sence of ionizing radiation. We explicitly include both reac- 
tions in Eq. ( |lO| ) via the term S''. An expression for this coeffi- 
cient can also be found in Arnaud & Rothenflug (1985 ). 



2.2.2. Collisional excitation of metastable states 

Neglecting the effects of radiative excitation and stimulated 
emission, we solve the equations of excitation equilibrium for 
the population densities Ny^ (sums over all values "j" for which 
the conditions under the summation signs are fulfilled): 



iVe E *j+ ^ ^kj 

E ^j9jk + E ^j^jk, (12) 

together with the condition A'j = A'gcs. We use the formu- 
lae for the activation and deactivation coefficients given in e.g. 



Osterbrock(1989): 



qi2 -- 
and 

921 = 



,.63.10-«^T-i/2 



exp 
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,63- IQ-^—T-^/'^, 

UJ2 



(13) 



(14) 



where £7x2 denotes the collision strength for the transition 1 
2, LUi and uj2 the statistical weights of both states involved and 
AEi2 the energy difference between them. For the ili2 we use 
the tables given in Osterbrock (1989). 



2.3. Balmer lines 

Our neglect of line absorption of Balmer photons by hydrogen 
is justified as long as the density of Ly^ photons is sufficiently 
low to insure that the hydrogen 2p state is not significantly pop- 
ulated. This is equivalent to the assumption that Ly^ photons 
generated in the nebula by recombination either are quickly 
destroyed, e.g. by dust absorption or by hydrogen Ly^ absorp- 
tion followed by 2-photon emission, or are able to escape suf- 
ficiently rapidly, e.g. by a random walk in frequency (Oster- 
brock 1961). The emission coefficient of the Balmer lines is 



given by: 



47r 



(15) 



The effective recombination coefficients used in this work 
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Fig. 2. Density, velocity and ionization structure of model A 
and C. Gray scale and black contour lines display the den- 
sity structure. These contour lines vary from logp = —13.0 
to logp = —19.5 in increments of Alogp = 0.5. The white 
contour lines mark the position of the ionization front and the 
arrows show the velocity field. The normalization is given at 
the upper right corner. 



2.4. Radiation from the central star 

As argued in Paper I the resulting UV spectrum of a star accret- 
ing material via an accretion disk is very uncertain. For sim- 
plicity we have assumed that the photospheric emission of the 
central source (star + transition zone) can be approximated by 
a black body of given temperature in the frequency range 
of interest (A < 100 nm). determines the "hardness" of the 
ionizing photons, thus affecting both the nebula temperature 
and the ionization fraction of oxygen and nitrogen. We use the 
same values for T* as in Papers I and II for the hydrodynamic 
models. 

Nevertheless, the successful spectral classification of the 
ionizing star i n the UCHll region G29.96-0.02 by Watson & 
Hanson ( 1997 ) gives rise to the hope that more information on 
the spectral properties of young, still accreting massive stars 
will be available in the future. 



were adopted from Hummer & Storey ( 1971 ) 



3. The Numerical Model 



O. Kessel, et al.: Photoevaporation of protostellar disks III. 



5 



G2 



Time [yr] 402+ 



v.,™ [Itm s"'] 400 



400 



200 



100 



0:32 



-1001/' 




G3 



-200 -100 100 200 300 

F [AJ] 



"Tims [yr] 4001 



v..™ [km =-'] 500 



4001 



300 



100 



G4 



-zoo -100 100 200 300 
R [AJ] 



Tima [yr] 3976 



v..m [Itm r'] 1000 



400 



100 



-20 



■ le" 



-14 



-12 




-20 



-16l^ 
a. 



■14 




-200 -100 100 200 300 
R [AJ] 

Fig. 3. Density, velocity and ionization structure of model G2, 
G3 and G4. Symbols and lines have the same meaning as in 
Fig. H except the black density contour lines, which are drawn 
down to log p = —21.5. 



3.1. Structure of the underlying models 

The underlying numerical models were calculated on five mul- 
tiply nested grids, each with 62 x 62 grid cells (see Yorke & 
Kaisig 1995 , Paper I, and Paper II). The spatial resolution of 
the finest grid was AR — AZ w 2 x 10^^ cm (R is the dis- 
tance to the symmetry axis, Z to the equatorial plane). Axial 
symmetry and mirror symmetry with respect to the equatorial 
plane were assumed for the models. The simulations were per- 
formed within a volume + Z^^^V^^ < 10^^ cm until a 
quasi-Steady state was reached. 

For the diagnostic radiation transfer calculations discussed 
here we use the final states of five simulations described in Pa- 
per II. Some of the relevant parameters of these simulations 
are given in Table |[ Figure ^ and Fig. ^ display the density 
and ionization structure as well as the velocity field of the se- 
lected models. Models A and C are the results of simulations 
with the same moderate stellar wind and the same radiation 
source. But in the simulation leading to model A the diffuse 
UV radiation field originating from scattering on dust grains 
was completely neglected. For that reason we got a higher pho- 
toevaporation rate Afph for model C. In Fig. ^ this is recogniz- 
able by the greater overall density in the ionized regions and by 
the higher velocity in the "shadow" regions of the disk in the 
case of model C. In order to investigate the variation of spec- 
tral characteristics with the stellar wind velocity we chose the 
models with the greatest wind velocities G2, G3 and G4. Fig- 
ure 1^ shows the increasing opening angle of the cone of freely 
expanding wind with increasing wind velocity. 



-12 3.2. Strategy of solution 

We use the model data to calculate the ionization structure and 
the level population. From the level populations we determine 
the emissivities of each line transition and the continuum emis- 
sion at each point within the volume of the hydrodynamic mo- 
del. For each viewing angle Q considered, we solve the time 
-20 independent equation of radiation transfer in a non-relativistic 
moving medium along a grid of lines of sight (LOS) through 
the domain, neglecting the effects of scattering: 



di, 
dr, 



(16) 



where the optical depth is defined as — J K^ds. Integra- 
tions were performed for a given set of frequencies, whereby 
the effects of Doppler shifts for the line emissivities were taken 
into account. The resulting intensities are used to determine 
SEDs, intensity maps and line profiles. Spectra are obtained 
from the spatial intensity distributions by integration, taking 
into account that each LOS has an associated "area". Depend- 
ing on the symmetry of the configurations could be utilized 
to minimize the computational effort (see Fig. Q). For the pole- 
on view (8 = 0°), for example, only a one dimensional LOS 
array need be considered. For the edge-on view (8 — 90°) 
lines of sight either through a single quadrant (continuum trans- 
fer) or through two quadrants (line transfer) are necessary. The 
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Fig. 5. Projection of a typical LOS trajectory (curved dashed 
line) onto the model data grid (solid lines). Temperature, den- 
sity, degree of ionization and velocity are defined at cell cen- 
ters. The small circles divide the LOS into subintervals; the 
source function Si, is evaluated at the location of the circles, 
chosen to lie on the intersections of the LOS with lines con- 
necting the grid cell centers. 

resolution of the central regions is enhanced by overlaying a 
finer LOS grid in accordance with the multiple nested grid 
strategy used in the hydrodynamic calculations. 

Each point in Fig. ^ corresponds to an LOS trajectory 
through the model. Mapping such a trajectory onto the (R,Z) 
model grid yields hyperbolic curves as displayed in Fig. ||. 
Beginning with a starting intensity {Iy{—co)), the solution of 
Eq. ( p^ is obtained by subdividing the LOS into finite inter- 
vals and analytically integrating over each interval assuming a 
sub-grid model (see below). 

3.3. Continuum radiation transfer 

If no discontinuities are present within the subinterval under 
consideration, we assume Si, varies linearly with r, i.e. 



Su{t) — si + (5* — ^Dly ; 

At 



(17) 



where i and i + 1 denote the starting and end points of the 
interval, respectively, and At is a mean optical depth over the 
interval: 



At = As. 



(18) 



With this formulation the sol ution o f Eq. (|l6|) over the entire 
interval is given by (see Yorke 1988): 



li+i = cxp(-At) + Si 



1 — exp(— At) 

a7 



exp(— At) 



For the cases considered here we choose /o = as the start- 
ing LOS intensity. For "proplyd"-type models (considered in a 
subsequent paper of this series) a non-negligible background 
intensity should be specified. 

3.4. Radiation transfer in emission lines 

For the transitions considered here the radiation field can be 
considered "diffuse" Ii, ^ and the contribution of spon- 
taneous emission dominates over line absorption and stimu- 
lated emission processes. After separating the source function 
S,y = Sc + 5'l and the intensity = Ic + h of Eq. (|l^ into 
the contributions of the continuum and the line, we obtain 



II = 



5'Lexp(-AT) 
V^Ai/th 



Ar 



exp 



v — v{t) 
Ai^th 



dT, (20) 



where cIt = kq As and 5l = cl/kc- Here v{t) is the Doppler- 
shifted frequency of the transition, Ai^th the Doppler width and 
the net source function integrated over the line. Assuming 
that v is linear in t over the whole interval yields the analytical 
solution to Eq. (pO|): 



, ^lAt 



(i>2 - I') At 
A^ 



[erf(r2)-erf(ri)],(21) 



where erf (y) = 2/ ^/tt exp{~t^)dt is the error function and 
Yi = {vi — v)l Ai^th a dimensionless frequency shift. 

The net source function S'l is calculated according to the 



algorithm suggested by Yorke (1988): 



~ erf(yM)-erf(yi) - ^x^ij^y-^d^~ 

— 77T7-\ 77T7tS\ H 7TT7\ n^r \ '^2, (22) 



erf(r2)-erf(ri; 



erf(r2)-erf(Yi; 



with 1m = (^1 + ^2)72. The line source functions S'l and 
^2 are calculated from the total line emission coefficient el as 
defined in Eq. ^ at the boundaries of the evaluated interval 
and from the continuum absorption coefficient: Si = eL.i/^c. 

If Ai> ^ Az^th, i e. there is negligible Doppler shift within 
the subinterval, the solution of Eq. (^0|) with vx ^ V2 ^ v vs, 
used: 



St 



V^Afth 



exp 



-X 2 
V — V 



At^th 



(1 - exp(-AT)) . (23) 



3.5. Treatment of ionization fronts 



The numerical models considered contain unresolved ioniza- 
tion fronts due to the coarseness of the hydrodynamic grid. At 
these positions jumps occur in the physical parameters and the 
solutions given by Eq. ^ and Eq. are poor approxi- 

mations. The exact location of the fronts within a grid cell are 
unknown; we assume they lie at the center of the corresponding 
interval. Our criterion for the presence of an ionization front is 




Fig. 4. Choice of lines of sight (LOS) and their associated areas for different viewing angles 9. Filled dots indicate the LOS used 
for the continuum calculations. Empty dots refer to the additional Lines of Sight necessary for the line profile calculations. 



a change in the degree of ionization Aa; > 0.1 between two 
evaluation points. 

For the continuum calculations Eq. (|l9|) is applied to each 
half interval with At = Kj Aa;/2. For the first half S = Si kept 
constant and for the second half S = Si+i is held constant. 
For the line calculations Eq. (^3|) is used with i> = z^i (1^2) and 

= Si {S2) for the first (second) half interval. 

3.6. Treatment of the central radiation source 

The central source is modeled by a black body radiator of tem- 
perature and radius R^. The integration along the line of 
sight through the center is started at the position of the source 
with the initial intensity 



(24) 



with A is the area associated with the central LOS. 
4. Results 

With the code described above we determined SEDs, contin- 
uum isophotal maps and line profiles for the forbidden lines 
[Nil] 6584, [On] 3726 and [Olll] 5007 as well as the Ha-Hne 
for the models introduced in Sect. 3J . Their relevant physical 
parameters are listed in Table ^ 

4.1. Continuum emission 

4.1.1. Spectral Energy Distributions 

Figure ^ shows the SEDs of model G2 presented in Paper II 
for different viewing angles <d. The spectra can be divided into 
three regimes dominated by different physical processes: 

1. In the frequency range from 10^ to 10^^ Hz the SED is 
dominated by the thermal free-free radiation in the ionized 
region around the dust torus. 



2. The iR-excess from 10^^ to 10^'' Hz is due to the optically 
thick dusty torus itself, which has a mean surface tempera- 
ture of about 250 K. 

3. Beyond 10^'' Hz the SED depends strongly on the viewing 
angle: if the star is obscured by the dusty torus then the 
free-free radiation of the Hll-region again dominates the 
spectrum, otherwise the stellar atmosphere shows up. Due 
to the uncertainties in the stellar spectra and the neglect of 
scattering by dust in Eq. (p^, which becomes more and 
more important with increasing frequency, a discussion of 
the SED beyond 10^'* Hz and comparison with observations 
are not useful. 



According to the analysis of Panagia & Felli ( |1975[ ), who 
calculated analytically the free-free emission of an isothermal, 
spherical, ionized wind, the flux density should obey a v^^^- 
law: 
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(25) 



Schmid-Burgk (1982) showed that this holds, modified by 
a geometry dependent factor of order unity, even for non- 
symmetrical point-source winds as long as p drops as R^^. Ad- 
ditionally he postulated that the flux density should hardly be 
dependent on the angle at which the object is viewed. In Fig, p 
we include for comparison the flux distribution given by Eq. 
for the photoevaporation rate Afph = 1.65 • 10^® M.Qyr~^ (see 
Paper II), T = 10 000 K and Wwind = 20 km s"^ derived from 



the line profiles in Sect. X.2. The slope of the SED in regime 



1 is slightly steeper in our results, because our volume of inte- 



gration is finite; Panagia & Felli (1975) derived their analyti- 
cal results by assuming an infinite integration volume. The flux 
is almost independent of 0, which is in good agreement with 
Schmid-Burgk (|1982|). The deviations between lO^'' and 10" 
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Table 2. Scattering coefficient ^I'Jfg* p ^ as well as parameters for the stellar wind (mass loss rate A/wind and velocity i^wind) and 
the ionizing source (stellar photon rate ^star and temperature Tcff) used in the calculations. The evaporation time scale icvap is 
calculated from tovap = M^isk/Mph with A/disk — 1-67 Mq. 



model 


^dust / P 
2 —1 

cm^g 


-^'-^wind 

lO^^Afoyr"^ 


t^wind 

km s^^ 


logjQ Sstar 

s-1 


Tcff 
K 


A/ph 
lO^^Afoyr^i 


^GVap 

lO'^yr 


A 
C 



200 


2 
2 


50 
50 


46.89 
46.89 


30000 
30000 


0.565 
1.35 


2.96 
1.24 


G2 
G3 
G4 


200 
200 
200 


2 
2 
2 


400 
600 
1000 


46.89 
46.89 
46.89 


30000 
30000 
30000 


1.65 
1.67 
1.64 


1.01 
1.00 
1.02 




fraquarey [H2] 
SED modal G4 



Fig. 6. Spectral Energy Distribution for model G2 and different 

e. 



Hz are due to the break in the R ^-law caused by the neutral 
dust torus. 

Figure ^ also includes the SED of a blackbody at T = 
240 K. The flux oc z^^-S in the far IR between 5 • 10" Hz and 
3 • 10^^ Hz is slightly steeper than the comparison blackbody 
spectrum, because the dust torus is not quite optically thick. 
With increasing Q the maximum shifts towards lower frequen- 
cies, because the warm dust on the inside of the torus is being 
concealed by the torus itself. 

We obtain qualitatively the same results for a number of 
models presented in Paper II. 

4.1.2. Isophotal maps 

We also calculated isophotal maps over the projected {S, T) 
grid of the sky for models C and G4 (Figs. The maps 
were convolved with a Gaussian point spread function (FWHM 
0. 3 for A = 6 cm, 0. 1 for A — 2 cm and 12 /im) in order 
to compare the numerical models with observations of limited 
resolution. The values (in percent) of the contour lines relative 
to the maximum flux per beam are 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 



20, 30, 40, 50, 60, 70, 80, 90 for A = 6 cm, and 2, 5, 15, 25, 35, 
45, 55, 65, 75, 85, 95 for A = 2 cm and 12/im. 

The difference in these two models is revealed most strik- 
ingly in the maps for A = 6 cm and 2 cm. The mass outflow 
of model C is more evenly distributed over its whole opening 
angle, whereas in model G4 most of the mass is transported 
outward in a cone between 8 = 30° and 70°. This leads to the 
X-shape of the corresponding radio maps for viewing angles 
> 60°. The high density in the region between star and disk 
for this model results in an optically thick torus in this region 
at A = 2 cm. Thus the contours for <d = 30° and 60° are not 
symmetric to the equatorial plane at T = . 

In the maps corresponding to 6 = 30°, A — 12/im there 
appears a peculiar horseshoe-like feature. It is generated by the 
hottest region of the dust torus, which is the innermost bound- 
ary with the smallest distance to the star. It can be seen as a ring 
in the maps for 6 = 0°. For 6 = 30° the part of the ring next 
to the observer is obscured by the dust torus, whereas the other 
parts are still visible. The beam width chosen by us allows the 
resolution of the ring and thus leads to the horseshoe-like fea- 
ture. For 9 = 60° only the most distant part of this hot ring is 
visible, leading to a maximum in the flux with a smaller spatial 
extent than for 8 = 30°. 

4.2. Line profiles 

Line profiles can provide useful information on the velocity 
structure in Hll regions. Because the thermal doppler broad- 
ening decreases with increasing atomic weight as A^^^"^, ions 
such as Nil and Olll are better suited for velocity diagnos- 
tics than H„. This can be seen in the line profiles we obtained 
(Figs. ^|l2|). They show the flux integrated over the whole area 
of the object including the disk, the evaporated flow and the 
cone of the stellar wind. In all cases the line broadening of 
several ten kms^^ is dominated by the velocity distribution 
of the escaping gas. For comparison, the rotational velocity of 
the dust torus Wjot — 13 kms~^ for the inner parts, and the 
thermal Doppler broadening from Eq. (§) at T = 10 000 K is 
wth — 9 kms^^ for Ha and ~ 2 kms^^ for Olll. 

Figures |l^ and |l^ show the line profiles for the models 
A (no scattering of H-ionizing photons) and C (calculated as- 
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Fig. 7. Isophotal maps of model C for different viewing angles and wavelengths as indicated. Assumed distance 300 pc. The 
circle in the lower right corner marks the FWHM of the point spread function. Values for contours see text. 
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Fig. 9. Lines [Nil] 6584 and Ha for models G2, G3 and G4, and different "viewing angles" 8. The models differ only in wind 
velocity Wwind- 



suming non-negligible UV scattering during the hydrodynam- 
ical evolution) from Paper II. UV scattering leads to stronger 
illumination of the neutral torus by ionizing radiation and thus 
to a higher photoevaporation rate in model C (by a factor of 
^ 2.3) compared to model A. Due to the higher density in the 
regions filled with photoevaporated gas, the lines for model C 
are generally more intense. In the case of the line [Oil] 3726 
one notices that the difference between the fluxes for different 
angles is the smallest of all transitions. Not only is the den- 
sity of the outflowing ionized gas higher for case C, but the 
charge exchange reactions discussed in Sect. 2.2.1 lead to the 
establishment of a non-negligible Oll-abundance even in the 
"shadow" regions not accessible to direct stellar illumination. 
These regions dominate the line spectra for all angles. 

Figures ^ and [ri| show the calculated line profiles for the 
models G2, G3 and G4, which only differ by the assumed stel- 
lar wind velocity Vvrmd, increasing from 400 kms^^ (G2) to 
1 000 km s~^ (G4). Comparing the results we find two features: 

1 . The intensity and overall structure of the line profiles con- 
sidered are almost independent of the stellar wind velocity 

'f^wind • 



No high- velocity component appears in the Unes due to the 
stellar wind. 



Both features can be explained by the fact that the den- 
sity of material contained in the stellar wind is much lower 
than the density in the photoionized outflow. Remembering that 
p = M / Anr'^v in a steady-state outflow and that the line emis- 
sivity El « p^, we can understand why the low expansion ve- 
locities of about 20 kms^^ (i.e. material close to the torus) 
dominate the spectra. The overall evaporation rates as well as 
the expansion velocities are almost equal for all three models, 
leading to very similar line profiles. It would be necessary to 
consider transitions which predominate in hot winds in order 
to detect this high velocity component and to find significant 
differences between these models. 

In spite of our assumption of optically thin line transfer, 
the profiles for Q ~ 90° are not symmetric. This is due to the 
dust extinction within the Hll-region. The receding material is 
on average further away from the observer than the approach- 
ing. The longer light paths result in stronger extinction of the 
redshifted components. 
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Fig. 10. Lines [Nil] 6584 and Ha for models A and C, and different "viewing angles" 9. In model A scattering of UV photons 
by dust was neglected, in model C included. Note the different scales on the abscissae. 



We are aware of the fact that the neglection of scattering 
by dust in Eq. (|l6|) may lead to serious errors in the calculated 
line profiles. We expect non-negligible contributions especially 
in the red-shifted parts of the lines due to light scattered by the 
dense, dusty surface of the neutral torus. This light was origi- 
nally emitted by gas receding from the torus. The resulting red- 
shift "seen" by the torus remains unchanged during the scatter- 
ing process and will thus lead to enhancement of the red-shifted 
wings of the lines. 

Nevertheless, we expect our qualitative results still to hold 
since the arguments mentioned above referring to the geometry 
of the underlying models are still applicable. 

5. Comparison with observations 

Although the cases presented here describe the situation of an 
intermediate mass ionizing source (8.4 Mq) with a circumstel- 
lar disk, many of the basic spectral features can be generalized. 
Thus, it is interesting to compare and contrast our results with 
observed UCHlls, even though many of the central sources are 
presumably much more massive. 

The collections of photometry data in the catalogues of 
Wood & Churchwell ( |1989| ) and Kurtz et al. ( |1994| ) show that 
the SEDs of almost all UCHll-regions possess roughly the 
same structure as the ones of the models discussed here: a flat 
spectrum in the radio and mm regime following a i/'^^-law due 
to free-free emission and an IR excess originating from heated 
dust exceeding the free-free emission by 3 — 4 orders of 
magnitude. A closer inspection shows that the dust tempera- 
tures are by an order of magnitude lower in most of the ob- 
served sources when compared to our models. This may be 
an indication that the disks are being photoionized by a close 



companion in a multiple system rather than the central source. 
Alternatively, the emitting dust could be distributed in a shell 
swept up by the expanding Hll-region and thus further away 
from the exciting star than for the cases discussed here. The 
large beam width of IRAS and the tendency of massive stars 
to form in clusters also make it likely that the IR fluxes be- 
long to dust emission caused by more than one heating star. 
Objects of this type would appear as "unresolved" in the maps 



presented by Wood & Churchwell (1989) for distances larger 
than ^ 300 pc. Due to the cooler dust, the "unresolved" ob- 
jects cannot be explained by the models of circumstellar disks 
around UV luminous sources specifically discussed in this pa- 
per. Certainly the cometary shaped UCHlls can be explained by 
a disk being evaporated by the ionizing radiation of an exter- 
nal star and interaction with its stellar wind. Numerical models 
dealing with this scenario will be presented in the next paper of 
this series. 



5.1. MWC349A 

A commonly accepted candidate for an evaporating disk 
around a young massive star is MWC 349 A. Its radio con- 
tinuum flux obeys the v^-^-\aw up to A = 30 cm (see collec- 



tion of photometry results in Thum & Martm-Pintado 1993). 
Radio maps obtained by Cohen et al. ( 1985| ) and White & 
Becker ( 1985 ) show an extended ionized bipolar outflow with a 
peculiar X-shape (Fig. pj| ). In the center Leinert (1986 ) finds an 
elongated, dense clump with T ^ 900 K, with optically thick 
IR emission, which makes MWC 349 A one of the brightest 
IRAS sources. The elongated structure shows an almost Kep- 
lerian velocity profile along its major axis, perpendicular to the 
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Fig. 11. Lines [Oil] 3726 and [Olll] 5007 for models G2, G3 and G4, and different "viewing angles" 6. The models differ only 
in wind velocity Wwind- 



outflow axis (Thum & Martm-Pintado 1993). This leads to the 
assumption of a neutral dust torus around t he cen tral star with a 
small outer radius < 100 AU. Kelly et al. ( 1994 ) estimated the 
extinction towards this early-type star to be Ay — 10.8. 

The SED of MWC 349 A shows all the features which 
we find for our models. The extinction and the geometry of 
the outflow, as well as the lack of a high-velocity compo- 
nent in the Une profiles (Hartmann et al. 



1980, Hamann & 



Simon 1988), could be explained by a model with fast stel- 
lar wind presented in this paper, assuming a viewing angle 
8 ~ 90°. On the other hand, the high dust temperatures in 
the torus, Td ^ 800 K, and the extremely high mass loss rate 
in the ou tflow, M = (1.16 ± 0.05) x 10"^ M© yr^^ (Cohen et 
al. 1985 ) remain puzzling and need clarification by a numerical 
model after the method described in this series but especially 
"tailored" to suit the needs of MWC 349 A. 



6. Conclusions 

In this paper we showed that the models of photoevaporat- 
ing disks around intermediate mass stars cannot explain the 
large number of "unresolved" UCHll's observed by Wood & 



Churchwell ( |1989| ) and Kurtz et al. ( |1994| ), because the inferred 
dust temperatures of these objects are in most cases an order of 
magnitude lower than those obtained in the numerical models. 
But the question remains whether the disks of more massive 
stars than considered here could be responsible for the high 
abundance of the "unresolved" UCHll's. Disks around close 
companions of massive stars should be treated in greater detail. 
If we assume that circumstellar disks are the rule in the process 
of star formation, the simplicity and straightforwardness of the 
model make it favorable compared to alternative suggestions. 
The extremely high radiation pressure in the vicinity of mas- 
sive stars could lead to a larger distance between star and disk 
and thus to smaller dust temperatures. 

Another important result of this work is that the profiles of 
forbidden lines in the optical for the models G2, G3 and G4 
with wind velocities of 400 — 1000 kms^^ are almost inde- 
pendent of Wwind- This is due to the fact that the mass loss rate 
and velocity of the evaporated disk material is not affected by 
its interaction with the stellar wind, but by the rate of ionizing 
photons and the peculiar shock structure which is very similar 
in the numerical models (see Fig. |[). Since the total mass loss 
rate is dominated by the evaporated component with low veloc- 
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Fig. 13. VLA-maps of MWC 349. Left: A = 6 cm (Cohen et al. |1985| ), FWHM 
90% of maximum flux 16,59 mJy beam^^. Right 
15, 25, ...95% of maximum flux 156 mJybeam^^ 
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ity, and the emission is proportional to p^, the intensity in the 
lines does not depend on details of the stellar wind and their 
profiles show the same low-velocity components. 

Treatment of optically thick line emission and scattering 
effects is not possible with the method presented above. In 
order to compare non-LTE-effects like masing lines, which 
are observed in various objects related with the formation of 
massive stars, one has to refer to different methods, e.g. the 
Monte-Carlo method presented by Juvela ( 1996| ). This would 
immensely help us in our understanding of the process of for- 
mation and evolution of massive stars. 
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